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QUASIINTERPOLANTS  AND  APPROXIMATION  POWER 
OF  MULTIVARIATE  SPLINES 


CARL  de  BOOR1 

Center  for  the  Mathematical  Sciences 
610  Walnut  Street 
Madison  WI  53705  USA 

ABSTRACT.  The  determination  of  the  approximation  power  of  spaces  of  multivariate  splines  with 
the  aid  of  quasiinterpolants  is  reviewed.  In  the  process,  a  streamlined  description  of  the  existing 
quasiinterpolant  theory  is  given. 


1.  Approximation  power  of  splines 

"  '  ■  •  p.  -•  '  , 

•f  >  j\  r*v  J  '  * 

'  I  begin  with  a  brief  review  of  the  approximation  power  of  univariate  splines  since  the 
techniques  for  its  investigation  are  also  those  with  which  people  have  tried  to  understand 
the  multivariate  setup.  (That  may  in  fact  be  the  reason  why  we  know  so  little  about  it.)  -L 
/-will  then  briefly  discuss 'three  examples  to  illustrate  some  basic  limitations  of  the  standard 
univariate  approach.  '  (^i'r  ^ _ . 

Let  S  :=  Sk, t  be  the  univariate  space  of  splines  of  order  k  with  knot  sequence  t. 
This  means  that 

S  :=  span  $ 


with 


$  :=  (v>i) 


n 

t=P 


<Pi  :=  . .  .,ti+k)  the  normalized  B-spline  for  the  knots  f,-, t,  and  t  := 

a  nondecreasing  real  sequence.  This  definition  of  a  spline  space  is  taylor-made  for  the 
consideration  of  its  approximation  power,  since  the  B-spline  basis 


n 

$  :  £oo(ra)  —  S  :  c  •-+  $c  :=  ifi  c(i) 

i=i 


is  so  well-behaved.  (I  have  found  it  convenient  to  identify  the  sequence  (<^j, . .  .,v?n)  with 
the  map  c  >-*  We  consider  specifically  approximation  from  S  to  X  :=  C([a,  i]), 

with 

[a, 6]  ;=  [t fc ?  tn+i] 
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the  interval  of  interest.  In  the  corresponding  norm 


ll/ll  •'=  {sup  |/(ar)|  :  a  <  x  <  h}, 


the  basis  (map)  $  satisfies 

||<I>||  =  sup  ||*c||/||c|U  =  1 

C 

(the  result  of  the  fact  that  the  </?,•  are  nonnegative  and  sum  to  1,  i.e.,  form  a  partition  of 
unity).  Since  any  linear  map  Q  on  X  into  S  has  the  form  Q  =:  for  suitable  linear 

functionals  A,,  it  follows  that 


||0||<||A||:=max||A,-||.  (1.1) 

The  (fi  have  local  support.  Therefore,  such  a  map  Q  is  local  to  the  extent  that  the  A i  are 
local.  In  [B68i],  I  chose 

supp  Aj  C  supp^t-  =  [t,-,ti+fc], 

and  will  take  this  to  be  the  meaning  of  the  statement  lQ  is  local’.  This  implies,  more 
precisely  than  (1.1),  that 


||Q/||[t.,t.+J]  <  ||A||  ||/||, 


+*  J 


(1.2) 


The  only  additional  feature  needed  to  make  Q  a  useful  approximation  scheme  from 
5  is  positive  (polynomial)  order  where  by  the  (polynomial)  order  of  any  such  map  Q  I 
mean  the  largest  integer  m  for  which 


<3  =  1  on  7r<m  :=  span(()J)j<m. 


(Strictly  speaking,  our  Q  is  only  defined  on  C([a,6]),  hence  I  should  talk  about  (7r<m)|[a  i)] 
instead  of  7r<m,  but  pedantry  can  be  overdone.) 

It  follows  that,  for  any  /  and  any  p  €  7T<m  and  any  /  := 

(/-<J/)(/)=((i-  OH/ -?))(/), 


hence 

ll(/  -  QW) II  <  111  -  <?||  ll(/  -  p)(ti+,_*,ti+*ilK 

with  |]1  -  Q||  boundable  by  1  4-  ||0||,  hence  by  1  +  ||A||.  Consequently, 


dist(/,5)  <  (1  +  !|g||)  rnaxdist(/,7r<m)[<i+J_titi+tJ.  (1.3) 

In  particular,  for  a  sufficiently  smooth  /, 


dist(/,T<m)(c,c+h,  =  0(i|Z?m/||hm), 


therefore 

dist(/,  5)  =  O(||0||  \\Dmf\\  |t|m),  (1.4) 
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with 


|tj  :=  max  A 


This  formulation  is  careless.  Offhand,  the  knot  sequence  t  enters  here  not  only  in  the 
meshsize  |t|  but  also  in  the  ||Q||,  since  the  actual  construction  of  suitable  Q  is  bound  to 
involve  the  knot  sequence.  There  are  many  ways  to  construct  local  Q  (choose,  e.g.,  each  A; 
to  be  evaluation  at  some  point  in  the  support  of  <?{).  There  are  also  many  ways  to  construct 
Q  of  any  order  <  k  (choose,  e.g.,  spline  interpolation,  or  least-squares  approximation  from 
S).  The  existence  of  a  Q  which  is  local  and  of  order  k  was  first  proved,  for  odd  k  and 
for  X  =  Ck([a,b])  and  without  the  aid  of  B-splines,  in  [Bi67],  with  a  corrected  version, 
which  also  covered  even  k,  to  be  found  in  [B682] -  For  the  present  situation,  i.e.,  for  X  = 
C([a,6j),  such  Q  were  first  constructed  in  [B681] .  The  Aj  there  were  even  chosen  to  be 
linear  combinations  of  point  evaluations.  Finally,  it  was  shown  there  that  the  A;  could  be 
so  chosen  that  ||A||  was  boundable  independent  of  t  (depending  only  on  the  order  k  of  the 
spline  space).  The  essential  part  of  the  argument  in  [B681]  is  the  observation  that  the 
B-spline  basis  is  locally  well  conditioned,  i.e.,  that 

1/dfc  :=  infmindist(v5,-,span(v7j)Jy<)tti,tl+k]  >  °- 

By  Hahn-Banach,  this  entitles  one  to  believe  in  the  existence  of  A*  on  A”  dual  to  $,  i.e., 
satisfying  A npj  =  with  suppA;  C  [£i, <,-+*,],  and  with  ||A||  <  d^.  The  resulting  Q  is 
therefore  not  only  of  order  k,  it  is  actually  the  identity  on  all  of  S,  i.e.,  it  is  a  linear 
projector  onto  S.  Since  7r<fc  C  5,  this  provides  the  bound 

dist(/,  5)  =  0(dfc||£)fc/||)|t|fc,  (1.5) 

in  which  the  order  term  0(- )  is  independent  of  the  knot  sequence  t. 

While  [B681]  contains  only  a  recipe  for  the  construction  of  suitable  A,,  [BF73]  estab¬ 
lishes  the  formula 

A,-  :  /  -  J2(-D)k-l-jMT,)(D]f)(Ti)  (1-6) 

i<k 

1>i  '■=  (<.+1  -  •)  •  •  •(<«+*- 1  -  0/(fc  -  1)! 

for  the  dual  functionals  for  the  normalized  B-splines  (with  r*  arbitrary).  While 

one  may  object  to  the  use  of  derivative  information  here,  the  formula  makes  it  easy  to 
extend  the  coordinate  functionals  of  the  B-spline  basis  to  linear  functionals  with  more 
desirable  characteristics  (such  as  applicability  to  Lp  functions  [B74],  or  employing  only 
function  values  [LS75]). 

On  the  other  hand,  one  may  choose  to  ignore  the  formula  (1.6)  and  choose  A,  explicitly 
so  that  Q  =  1  on  x<m.  This  requires  each  A;  to  be  an  extension  from  7r<m  only,  of  the 
linear  functional  which  associates  with  each  p  €  7r<m  its  ith  B-spline  coefficient.  If  each 
extension  happens  to  have  its  support  in  an  interval  containing  no  knots  in  its  interior  and 
if  m  =  fc,  then  the  resulting  Q  will  be  the  identity  on  all  of  S. 

The  specifics  just  discussed  provide  an  instance  of  the  following  more  abstract  situation: 
We  wish  to  approximate  from  a  given  directed  family  ( Sh )  of  function  spaces  in  C(G),  for 
some  compact  domain  G  C  IRd.  Specifically,  we  are  interested  in 

dist(/,  Sh) 
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as  a  function  of  the  parameter  h  whose  corresponding  meshsize  \h\  we  think  of  as  going 
to  zero  (and  it  is  in  this  sense  that  we  think  of  (Sh)  as  ‘directed’).  We  are  able  to  exhibit  a 
good  quasiinterpolant  Qh  of  polynomial  order  m,  i.e.,  a  linear  map  Qh  on  C(G)  into 
Sh  of  (polynomial)  order  m  (meaning  that  Qh  =  1  on  7r<m)  which  is  uniformly  local  in 
the  sense  that 

I (Qhf)(x)\  <  const||/|Br|M(l)|| 

for  all  x  €  G,  for  some  constants  const  and  r,  and  with  Br(x)  the  ball  of  radius  r  around 
x.  This  implies  the  error  bound 


dist (f,Sh)  <  const/|h|m 

valid  for  all  ‘sufficiently  smooth’  /  and  so  establishes  the  approximation  order  from  (Sh) 
to  be  (at  least)  m.  Further,  we  are  able  to  establish  m  to  be  the  (exact)  approximation 
order  from  (Sh)  by  exhibiting  a  particular  ‘sufficiently  smooth’  function  /  for  which 

dist(/,Sh)  ^  o(|/i|m). 

Finally,  we  are  able  to  choose  Qh  to  be  even  a  projector  onto  Sh- 

This  abstraction  has  motivated  much  of  the  work  on  the  approximation  power  of  mul¬ 
tivariate  splines.  But,  before  starting  that  discussion  here,  I  want  to  bring  three  cautionary 
examples. 

Example  1  This  example  appears  in  [DR8x].  It  concerns  the  space  Sh  spanned  by 
the  /iZZ-translates  of  the  piecewise  linear  function 

fx+1,  0  <  x  <  h\ 

^  •  X  ^  {  0,  otherwise. 

Thus  Sh  consists  of  certain  piecewise  linear  functions,  with  breakpoint  sequence  h7Z,  but 
the  only  constant  function  it  contains  is  the  constant  0.  In  particular,  it  is  not  possible  to 
construct  a  quasiinterpolant  of  positive  order  for  it.  Nevertheless,  the  approximation 

Qnf.=  Y,  VH(--i)fU) 

iei’iz 

has  the  error 

j£hTL  jZhTL 

with  Xh  the  characteristic  function  of  the  interval  [0,  /i[,  hence  ||y*  -  <ph\\  —  h.  Therefore 

11/  -  Qhf\\  <  w/(h)  +  !l/l|h, 

with  Uf  the  modulus  of  continuity  of  /.  It  follows  that  Qhf  converges  to  /  uniformly  in 
case  /  is  uniformly  continuous  and  bounded. 

This  example  strongly  stresses  that,  in  the  earlier  argument  involving  good  quasiinter- 
polants,  positive  order  is  too  strong  a  condition.  In  fact,  as  is  evident  from  the  argument’s 
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details,  it  is  sufficient  to  demand  that  the  good  quasiinterpolant  Qh  be  of  positive  local 
order  m,  meaning  that 

Qhf  =  /  +  0(11/11  \h\m) 

on  any  ball  of  radius  0(|/i|)  and  for  any  /  6  ff<m-  A  sufficient  condition  for  this  is  that 
Qh  =  1  on  a  /^-invariant  space  Fh  of  entire  functions  which  is  locally  close  to  7r<m  in  the 
sense  that  its  ‘limit  at  the  origin’,  (Fh)i,  contains  7r<m.  (Here,  F;  :=  span{/|  :  /  e  F}, 

with  fi  the  first  nonzero  term  in  the  expansion  /  =  /o  +/i  H - of  /  into  homogeneous 

polynomials  fj  of  degree  j,  all  j.) 

This  example  also  illustrates  the  limits  of  the  Strang-Fix  conditions  (see  Section  4).  For 
it  shows  that  (Sh)  has  positive  approximation  order  even  though  none  of  the  Sh  contains 
7r0.  In  fact,  =  {0}. 

Example  2  This  disturbing  example  comes  from  [BH83].  The  space  S  :=  of 
C1-cubics  on  the  three-direction  mesh  A  contains  all  cubic  polynomials.  It  even  contains 
them  locally  in  the  sense  that  any  cubic  polynomial  on  one  of  the  triangles  of  the  partition 
A  can  be  extended  to  an  element  of  S  with  compact  support.  However,  the  expectation 
raised  by  this  that  the  approximation  order  from  ( Sh  •'=  ,,A)  is  4  will  be  disappointed. 

There  are  polynomials  in  7 r4  for  which  dist(/,  ir\  ^  o(h3).  This  indicates  that  it  is  not 
sufficient  to  find  out  which  polynomials  are  contained  in  Sh-  One  needs  to  know  that  they 
are  contained  in  5*  in  a  local  and  stable  way. 

Example  3  In  [BHS87],  Hermite  interpolation  to  planar  curves  by  parametric  piece- 
wise  cubic  curves  is  investigated.  The  curves  being  piecewise  cubic,  one  would  expect  an 
approximation  order  of  4,  i.e.,  an  error  of  0(|/i|4),  with  \h\  a  measure  of  the  spacing  of  the 
length  of  the  cubic  pieces.  But,  in  fact,  the  scheme  described  there  is  shown  to  approximate 
smooth  curves  (without  inflection  points)  to  0(|fi|6).  Further,  an  example  is  given  there 
of  a  convex  smooth  curve  with  a  flat  spot  to  which  the  scheme  approximates  only  within 
0(|/i|4).  Finally,  [S8x]  shows  that  the  scheme  can  be  appropriately  modified  to  approximate 
to  smooth  curves  with  simple  inflection  points  to  within  0(|/i|6). 

Now  add  to  this  the  fact  that  the  approximation  scheme  used  is  nonlinear.  This 
means  that  the  basic  trick  of  the  above  argument,  viz.  the  introduction  of  a  local  polynomial 
approximation  to  /,  is  not  readily  applicable.  Still,  the  only  ready  means  for  estimating 
approximation  order  is  the  introduction  of  a  local  polynomial  approximation.  Now  add  to 
this  the  additional  difficulty  that  there  is  no  natural  way  of  measuring  the  distance  between 
curves,  other  than  their  Hausdorff  distance  which  is  apt  to  set  up  a  not  very  smooth  map 
between  the  points  of  the  two  curves. 

The  essential  facts  in  the  proof  of  the  0(|/i|6)  error  turned  out  to  be  that  (i)  the  scheme 
is  local,  and  (ii)  each  cubic  curve  piece  matches  the  given  curve  in  six  independent  pieces  of 
information.  This  is  sufficient  to  show  that,  for  sufficiently  small  |/i|,  the  difference  between 
each  cubic  piece  and  the  given  curve,  measured  in  a  suitable  local  direction  (i.e.,  interpreting 
both  as  (graphs  of)  functions  in  some  suitable  local  coordinate  system),  is  0(|/i|6),  with  a 
constant  that  can  be  bounded  in  terms  of  the  local  radius  of  curvature  of  the  given  curve. 

All  in  all,  it  is  again  the  use  of  a  local  and  stable  scheme  of  some  positive  order  that 
supplies  the  approximation  order,  but  the  argument  is  much  harder,  and  is  unsatisfactory 
precisely  because  it  takes  recourse  to  functions.  It  seems  pretty  hopeless  to  try  to  settle 
approximations  of  surfaces  by  surfaces  in  this  way.  Hence,  although  this  is  a  meeting 
on  computing  with  curves  and  surfaces,  I  will  stick  to  the  approximation  of  multivariate 
functions. 
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2.  Quasiinterpolant  construction  for  a  scale 

We  now  commence  the  discussion  of  quasiinterpolants  and  approximation  order  for  multi¬ 
variate  spline  spaces. 

The  simplest  model  for  a  family  (5\)  of  approximating  spaces  is  that  of  a  scale,  i.e., 

Sh.  '■=  (?h{S) 

for  some  fixed  space  5,  with 

Oh.f  :  x  >—  f(x/h) 

and  h  >  0.  The  simplest  nontrivial  model  for  S  investigated  is  that  of  the  span  of  the 
integer  translates  of  a  compactly  supported  function.  This  means  that  S  is  taken  to  be  of 
the  form 

S(<p)  :=  ran  tp*  :=  {<p  *  c  :  c  €  C2  }. 

This  description  makes  use  of  the  convolution 

c  :=  ¥>(•  -  a)c(a) 

a€  TL* 

of  the  compactly  supported  function  <p  with  the  complex-valued  sequence  or  meshfunction 
c,  i.e.,  c  :  7Zd  — <■  C.  Since  tp  is  compactly  supported,  the  infinite  sums  ip*c  converge  trivially 
uniformly  on  compact  sets,  and  it  is  in  this  sense  that  I  mean  to  interpret  them. 

I  now  consider  the  problem  of  constructing  a  good  quasiinterpolant  of  polynomial  order 
m  for  S(<p).  Following  [BH82],  it  has  become  customary  to  construct  such  a  quasiinterpolant 
in  the  form 

<3  (2-1) 

a 

for  some  suitable  local  bounded  linear  functional  A.  With  the  notation 

V  *'  f  ■=  V  *  f\  =  52  v(-  ~  “)/(«*) 

a 

for  the  semi-discrete  convolution,  which  uses  the  helpful  abbreviation 

f\  :=  f]z*> 

and  the  notation 

A  fix*-*  A /(•  +  i), 

we  can  write  (2.1)  more  simply  as 

Q  :  f  ip*'  A/. 

With  both  A  and  tp  compactly  supported,  the  corresponding  family  Q/,  :=  <ThQtJ\/h  is 
uniformly  local,  hence  establishes  approximation  order  m  for  the  scale  (S(<p)/i)  provided 
Q  =  1  on  7T<m. 


6 


r 


Assuming  that  7r<m  C  5(y?),  this  raises  the  question  of  just  how  one  might  construct 
Q  of  the  form  (2.1)  so  that  it  is  of  polynomial  order  m.  Consider  the  following  slightly  more 
general  problem:  Construct  A  so  that  Q  of  (2.1)  is  the  identity  on  the  finite-dimensional 
linear  space  F. 

Such  a  space  F  must  necessarily  lie  in  5(y>).  Further,  since  Q  is  to  be  the  identity  on 
it,  we  might  as  well  assume  that  F  is  E-invariant,  meaning  that 

EaF  C  F  V(q  €  7Ld), 


with  E  the  shift,  i.e., 

r/:iH/(i  +  a). 

For,  the  shift  commutes  with  <p*'  and  with  A,  hence  with  Q,  and  this  implies  that  Q  =  1  on 
the  smallest  E-invariant  space  containing  F.  (This  would  imply  that  F  is  even  F-invariant, 
i.e.,  invariant  under  differentation,  in  case  F  C  tt;  cf.  [B87]). 

It  follows  that  F  is  also  -invariant,  as  is  any  E-invariant  subspace  of  S(<p),  since 

=  V(/ =  *>*<:  €$(*>)).  (2.2) 

(Here  are  the  details  from  [B87] :  <p  *'  f  =  <p  *  (p>\  *  c)  —  ip  *  (c  *  y?|)  =  (y>  *  c)  *  qp|  =  /  *'  ip). 
This  means  that  the  linear  map 

T  :=  p  *'  j f 

carries  F  into  itself. 

Assume,  in  addition,  that  T  is  1-1,  hence  invertible.  Then  any  A  which  is  an  extension 
of  the  linear  functional 

A0  :  F  -  C  :  /  ~  (T~1/)(0)  (2.3) 

provides  Q  =  1  on  F,  since  for  such  a  A  and  for  any  /  £  F  and  any  a  £  7Zd , 

A/(-  +  «)  =  (T-lF°/)(0)  =  (E°T-1/)(0)  =  (T-1/)(a) 

(using  the  fact  that  T  commutes  with  E,  hence  so  does  T-1),  and  therefore  A/  =  T-1/  on 
7Zd,  consequently  Qf  =  TT~X  f  =  /. 

If  the  translates  of  p>  are  linearly  independent,  then  every  /  €  F  has  a  unique  repre¬ 
sentation  in  the  form  y?  *  c,  hence  Q  =  1  on  F  only  if  A  is  an  extension  of  the  functional 
A0  =  [0]T-1.  Here  and  below,  [x]  denotes  the  linear  functional  [x]  :  /  *—  f(x)  of  point- 
evaluation. 

Since  the  action  of  Q  on  F  is  decided  entirely  by  the  values  the  linear  functionals 
/  ►—  A/(-  +  a)  take  on  F,  it  is  easy  to  enlarge  the  class  of  available  quasiinterpolants  by 
considering  (as  already  cited  earlier  from  the  univariate  literature,  but  see  also  [CD88]  and 
[BeR8x]),  more  generally, 

Q  :  /  ~  “  Q)  A<*/ 

a 

with  Aa  =  [0 }T~lEa  =  [a]T-1  on  F. 

This  leaves  the  question  of  just  how  one  might  actually  construct  suitable  A.  I  take 
up  this  question  in  Section  5,  after  a  discussion  of  the  Strang-Fix  conditions.  But  it  is 
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convenient  here  to  discuss  quickly  the  special  case  when  F  is  a  polynomial  space.  Since,  by 
(2.2),  for  any  p  €  n  ft  5(<p), 


<p*'p  =  p*'<p€p^  <f(a)  +  7r<degpi 

Q 


T  is  1-1  if  and  only  if  <p(a)  /  0.  With  this  assumption,  we  may  assume  without  further 
loss  of  generality  that  <p  is,  in  fact,  normalized,  i.e.,  that 

=  L  (2.4) 


Then,  the  linear  functional 

r  :=  [0]v?*'  =  ^  -<*)[<*] 

a 

is  a  finite  weighted  sum  of  function  values  and  takes  the  value  1  at  the  constant  function 
1.  This  makes  it  possible  to  construct  numerically  its  associated  Appell  sequence  ( pa ) 
(see  the  appendix  for  complete  relevant  details  on  Appell  polynomials).  This  sequence  is 
characterized  by  the  fact  that  V(a,/3)  rDapl 3  =  This  implies  that  pa  €  na,  hence  the 

sequence  can  be  constructed  numerically  from  the  numbers  <p|  by  recursion: 

pa  =  n*-5>ur~v*-  (2.5) 

0<a 

Here,  [|°  :  x  ►—  xa /a\  is  the  normalized  power  function.  The  Appell  sequence  is  of 
interest  here  since  ([B87]) 

v>*'pa  =  ir 

for  all  [f0,  €  5(<p)-  This  implies  that,  for  all  such  a, 

Ao0°  =  Pa(  0), 

thus  providing  us  with  the  matrix  representation  of  Ao  with  respect  to  the  normalized  power 
basis  (if  any)  for  F.  In  particular,  since  [0]L>a[J>i  =  <50_p,  this  provides  immediately  the 
extension 

A  :r_  53Pj(0)(0]D° 

or 

of  A0. 

If  F  does  not  have  a  normalized  power  basis,  or  is  not  known  exactly,  then  it  is  useful 
to  observe  that  T  is  the  restriction  of  <p|*  to  F,  where  we  now  interpret  <pj  as  the  linear 
functional  f  >-*  =  r /(-•),  hence  (v?| */)(*)  =  <P| /(*-•)  =  £0  f(x -a)<r(a) 

and 

V\*Po  =  Ua  V(q  €  7Ld). 

The  assumption  (2.4)  guarantees  that  is  invertible  on  any  F-invariant  polynomial  space 
and,  in  particular,  on  any  superspace  Xfc  for  F.  This  implies  that  ^|a|<fc P<*(0)  [0]Frt 
extends  [0](v?|*)|^1 ,  hence  also  extends  A0. 
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This  observation  also  makes  apparent  the  following  comment  concerning  the  numbers 
pa(0).  If,  as  is  certainly  the  case  here,  r  is  a  compactly  supported  distribution,  then  (cf. 
Appendix) 

pc  =  i  -ii>ni/*,)(o) = "£  D“~’  [on  -  wra/ir).  as) 

Tf<Or 

This  means  that 

Pa(o)  =  [-tDr(i/^)(o), 

where,  in  our  case, 

=  <f\  =<P=Yl  iP(a)e~,0,i) 

a 

is  a  trigonometric  polynomial,  the  symbol  of  <p.  In  particular  cases,  it  might  be  easier  to 
compute  the  first  few  terms  in  the  Taylor  expansion  of  its  reciprocal  directly  rather  than 
by  the  recurrence  (2.5). 
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3.  Characterization  of  local  approximation  order 

The  Strang-Fix  conditions  originally  served  (see  [SF73]  and  references  there)  to  characterize 
the  ‘controlled’  approximation  order  from  the  scale  (S(<p)h).  The  surprising  aspect  of  this 
result  was  the  claim  that  having  7r<m  C  S(<p)  was  necessary  if  (S(p)k)  was  to  have 
approximation  order  m.  The  conclusion  drawn  from  this  (perhaps  too  eagerly)  was  that 
we  might  as  well  restrict  ourselves  to  piecewise  polynomial  <p.  Yet,  as  the  example  from 
[DR8x]  quoted  earlier  (and,  in  fact,  the  results  in  [DR8x]  on  the  approximation  power  of 
S(p)  when  <p  is  a  piecewise  exponential)  indicate,  there  are  perfectly  reasonable  ‘directed’ 
families  (Sh)  with  positive  approximation  power  which  contain  not  a  single  (nontrivial) 
polynomial. 

Following  a  suggestion  by  Babuska,  [SF73]  considered  the  more  general  problem  when 

5  =  S($)  :=  ^  S(v>) 

for  some  finite  set  $  of  compactly  supported  functions.  The  basic,  and  perhaps  surprising, 
result  is  that,  in  effect,  nothing  new  happens:  The  ‘local’  approximation  order  (defined 
after  the  statement  of  the  Theorem  below)  of  such  a  scale  is  the  best  that  can  be  had  from 
any  S(rp)  with  i p  in 

5X$)loc  :=  {]>^<p*cv  :  #suppcv  <  oo}. 
v> 

In  other  words,  even  for  5(4>),  the  ‘local’  approximation  order  is  realizable  by  a  quasiinter- 
polant  of  the  simple  form  ]T)0  E~a‘ipXEa  discussed  in  the  previous  section.  But,  in  contrast 
to  the  case  when  $  is  a  singleton,  there  is  at  present  no  computational  procedure  for  the 
construction  of  a  suitable  ip  or  for  the  determination  of  the  approximation  order. 

The  univariate  case  and  the  experience  with  the  simple  scale  (S(<p)/,)  related  in  the 
previous  section  gave  rise  to  the  hope  that  the  approximation  order  of  (5($)/,)  would  be 
the  largest  m  for  which  7r<m  C  5(4>).  Even  the  determination  of  such  an  m  would  be 
nontrivial,  but  less  involved  than  finding  a  best-possible  ip.  Unfortunately,  any  such  hope 
was  dashed  in  [BH83]  where  it  is  shown  that  the  approximation  order  of  the  scale  obtained 
from  the  space  of  C^-cubics  on  the  three-direction  mesh  is  only  3  even  though  its  subspace 
5($)  with  $  consisting  of  the  two  box-splines  M221  and  M 122  contains  ^<4. 

[SF73]  as  well  as  [DM84]  speak  of  ‘controlled’  Lp-approximation  order  and  mean  by 
that  the  largest  m  so  that,  for  all  sufficiently  smooth  functions  /, 

dist  „(/,$(*)')  <  const||Dm/||p(5r,(G))hm, 

with 

$(*){  •=  {<Th  £  <p  *  c„  :  £  ||cjp  <  const\\f\\p(Brfl(G))/hd/r>}  (3.1) 

for  some  const  and  r  independent  of  h  or  /.  Here,  the  distance  is  to  be  measured  in  the 
p-norm  on  the  underlying  domain  G.  Jia  gave  an  example  in  [J84]  to  show  that,  contrary 
to  the  assertions  in  [SF73J,  the  ‘controlled’  approximation  order  defined  this  way  for  the 
case  G  =  IR^  cannot  be  characterized  by  the  Strang-Fix  conditions.  Jia’s  example  does 
not  contradict  the  claims  made  in  [DM84]  (in  reliance  on  [SF73])  since  there  the  above 
‘controlled’  approximation  order  condition  is  assumed  to  hold  for  all  domains  G ,  with  const 
and  r  also  independent  of  G.  In  fact,  these  claims  are  verified  in  [BJ85]  where  the  following 
corrected  version  of  the  Strang-Fix  theorem  is  proved. 
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(3.2)Theorem.  Let  $  be  a  finite  collection  of  compactly  supported  essentially  bounded 
functions  on  IR^.  Then  the  following  conditions  are  equivalent. 

(i)  For  some  sequence  (ipa) |a|<m  in  span 

(ia)  ^o(O)  =  1,  xj)o  =  0  on  2n7Ld\0; 

(ib)  5Z/3<a  [  ~  =  0  on  2x2^0  for  0  <  |a|  <  m. 

(ii)  For  some  sequence  (^a)|a|<m  in  span$, 

[T  -  ^2  *'  [jf  e  7r<|a|  V(|a|  <  m).  (3.3) 

0<a 

(iii)  For  some  rp  6  5($)loc, 

[f -</>*' If  <=  7T<|ori  V(|a|  <  m).  (3.4) 

(iv)  For  all  p  £  [1,  col,  $  provides  local  Lp -approximation  order  m. 

(v)  For  some  p  €  [1,  oo],  $  provides  local  Lp -approximation  order  m. 

Several  comments  are  in  order. 

The  statement  provides  local  Lp-approximation  order  m’  is  meant  as  an  abbreviation 
for  the  condition  that,  for  all  /  6  L™(IRd), 

dist  p(/,  5($){)  <  const||£)m/||p  hm 

with 

sml  -  £  <p  *  c^  :  dist(j7i,  supp /)  >  r  =>  c^{j)  =  0}  (3.5) 

yjglt 

and  for  some  const  and  r  independent  of  h  or  /.  Here,  all  norms  are  the  Lp-norm  on 
JRd.  Note  that  the  claimed  equivalence  between  (iv)  and  (v)  makes  it  possible  to  drop  the 
qualifier  Lp-  and  talk  simply  of  the  local  approximation  order  provided  by  4>. 

The  proof  of  (i)=>(ii)  uses  Poisson’s  summation  formula,  in  the  manner  detailed  in  the 
next  section,  the  proof  of  (ii)=>(iii)  obtains  ^  in  the  form  )C|/3|<m  *  CP  w't^1  CP  finitely 

supported  sequences  for  which  fj0'  *  cp  =  0‘3  for  all  (3  <  7  :=  (m, . . . ,  m),  while  the  proof  of 
(iii)=>( iv)  uses  the  material  detailed  in  the  preceding  section.  The  main  point  of  [BJ85]  is 
a  proof  of  the  implication  (v)=>(i)  which  had  been  missing  in  the  earlier  literature.  The 
argument  uses  a  smooth,  compactly  supported  function  /  whose  Fourier  transform  satisfies 

7(0)  =  1,  Daf(x/h)  =  0(hm )  V(x  €  IRd\0,  |a|  <  m).  (3.6) 

(The  specific  function  /  used  happens  to  be  the  d-fold  tensor  product  of  the  centered  B- 
spline  of  order  m+  1,  hence  its  Fourier  transform  is  the  d-fold  tensor  product  of  Whittaker's 
sine  function  t  >-»  sin(l/2 V( ' / 2 ) ,  but  such  detail  doesn't  matter.)  Condition  (v)  provides 
an  approximation  fh  =  'rk  £)v9*c,p  to  /  which  has  compact  support,  uniformly  in  h.  and 
for  which 

:=/-A  =  0(/im) 

in  whatever  p-norm  we  happen  to  use.  It  follows  that  its  Fourier- Laplace  transform  ? :  :  >— 
fire*  e~'zz t(x)dx  satisfies 

|c(r)|  =  0(/im)  uniformly  on  1 1 1 rn  r||  <  const 
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for  some  const,  hence  Cauchy’s  formula  supplies  the  estimate 


||Ua?l|0o(IR‘i)  <  const  hm.  (3.7) 

Consequently,  the  function 

h=:J2v(h-K(h-) 

satisfies  the  conditions  (3.6)  (with  /  replaced  by  fh)  as  h  -*  0.  Here, 

vhv:z^hdY^e-izjchv(j) 

i 

is  2?r-periodic,  therefore 

[  -  wrw^j/h)  =  fcW  Y,  E  [  -  -  iD]a-0v^(O). 

V  0 

Combining  (3.6)  and  (3.7),  we  conclude  that 

lim  V(3(0)v»  =  l  (3.8) 

h— -*0 

and,  for  all  j  €  7Ld\ 0  and  all  |a|  <  m, 

Jim  E  E  l  ~  vW)I  -  iD\a-0v*{ 0)  =  0.  (3.9) 

<fi  p 

Define  now 

0-v  :=  E^-’ 

with  w  taken  from  the  orthogonal  complement  Wj_  of  the  space  W  of  all  ‘sequences' 
w  =  (inv,7)  for  which 

Jim  E  E  “Wl  “  iDVvv(0)  =  0. 

V  |7|<m 

We  may  choose  w  so  that  0o(O)  =  1,  since  0o(O)  =  0  for  all  in  _L  W  would  imply  that  the 
sequence  )>-*  <t>{0)8y  is  in  Vl7±±  =  W,  thus  lim/,_o  ^(0)^(0)  =  0  and  this  would 
contradict  (3.8).  It  follows  from  (3.9)  that  the  sequence  (^>,7)  1—  |  —  t£>]c,-"l'<JS(2:rj)  is  in 
W  for  all  |a|  <  m  and  all  j  €  27rZd\0,  hence 

E  l  -  iDfi>a-p{2TTj)  =  E  l  -  iDf  =  °’ 

0  0  V 

and  this  finishes  the  proof. 

The  proof  (i)=>  (ii)  =>(iii)  supports  the  stronger  claims  that  the  displayed  functions 
in  (ii)  and  (iii)  are  necessarily  in  x<a. 
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The  condition  actually  used  in  [SF73]  instead  of  (ii)  is  the  existence  of  some  sequence 
(V’o)iorl <rn  in  span  $  for  which 

[]"  =  £  •'  []'’  V(|a|  <  m).  (3.10) 

p<a 

Although  seemingly  stronger,  this  is  actually  equivalent  to  (ii),  as  can  be  seen  directly  as 
follows:  There  is  a  unique  linear  map  F  on  ir<m  which  carries  []a  to  Yjp<a  *'  Q  • 
By  (ii),  F  carries  v<m  into  itself  and  is  degree-preserving,  hence  must  be  1-1,  hence  must 
be  invertible.  This  is  all  happening  on  a  finite-dimensional  linear  space,  therefore  p(F)  =  1 
for  some  polynomial  p.  But,  for  any  polynomial  p ,  p{F)  carries  []a  to  Yjp<aXo,-p  *'  [] 
for  certain  x-y  expressible  as  linear  combinations  of  the  ifry. 

The  theorem  characterizes  the  local  approximation  order  from  the  scale  (5($)/!)).  The 
only  results  available  so  far  on  the  approximation  order  itself  have  been  obtained  by  produc¬ 
ing  an  upper  bound  on  the  approximation  order  which  happened  to  coincide  with  the  local 
approximation  order.  This  is  so  when  $  consists  of  a  single  box  spline  (cf.  [BH82]),  and  is 
often  so  when  $  consists  of  all  the  different  box  splines  in  with  A  the  three-direction 
mesh  (cf.  [J86]). 
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4.  The  Strang-Fix  conditions 


The  Strang-Fix  conditions  arose  in  the  characterization  of  the  approximation  power  of  the 
scale  ( Sh )  when  5  =  S(^>).  According  to  [SF73;  Theorem  I]  and  roughly  speaking,  the  scale 
(■SX^O/i)  has  approximation  order  m  if  and  only  if  ir<m  C  S(v>).  In  the  proof,  the  basic  tool 
is  the  Fourier  transform 

<p  :  £  Je-^(x)dX, 

with  fx  :=  (f,x)  =  YLj  the  scalar  product.  Since  < p  has  compact  support,  ip  is  an 

entire  function.  As  shown  in  [SF73],  the  following  conditions 

<£(0)  =  I  (4.1a) 

V(/3  <  a)  D0Cp  =  0  on  2xZd\0,  (4.1b) 

called  the  Strang-Fix  conditions  of  index  a  these  days  (cf.,  e.g.,  [C88]),  imply  that 

ira  :=  span  ( ()0)p<a  C  S(<p). 

[DM83]  prove  that,  more  generally,  any  affinely  invariant  subspace  P  of 

{p  6  Jr  :  p(D)(p  =  0  on  27r5Zd\0}  (4.2) 

is  necessarily  in  S(<p).  This  is  indeed  a  generalization  since  the  polynomial  space  ira  is,  in 
particular,  affinely  invariant,  i.e.,  scale-  and  translation-invariant.  Their  proof,  as  does  the 
argument  in  [SF73],  involves  the  semidiscrete  convolution 

9  *'/:=¥>  *(/,)=  £ 
pen* 

in  which 

f\  '■=  f\TZd- 

It  is  observed  in  [B87]  that  the  assumption  of  affine  invariance  can  be  weakened  to 
.E-invariance,  i.e.,  to  the  assumption  that 

V(a  €  ZZd)  E°P  C  P, 

with  E  the  shift,  i.e. 

E°f  :  x  >-*•  /(x  +  a). 

Explicitly,  [B87;  Prop. 2.2]  proves  that,  more  generally  (and  without  the  assumption  that 

*5(0)  *  0), 

{p  e  JT  :  (p*'p  =  p  *' <p]  (4.3) 

is  the  largest  E-invariant  subspace  of 

II*  :=  {p  €  7r  :  p(-iD)ip  =  0  on  27rZd\0}. 

Note  that  the  possible  lack  of  scale-invariance  (cf.  [R8x]  for  a  <p  for  which  (4.3)  fails  to  be 
scale-invariant)  forces  the  switch  from  p{D)  in  (4.2)  to  p(—iD)  here.  Note  further  that  (4.3) 
is  the  collection  of  polynomials  in  the  kernel  of  the  commutator  of  <p,  i.e.,  of  the  map 

/  [<p\f]  ••=  <P  *'/“/*'  V5 

singled  out  in  [CJW87]  as  an  object  of  interest. 

Finally,  [BR8x]  prove  the  following  more  general  and  suggestive  result,  in  which 

is,  by  definition,  the  largest  E-invariant  subspace  of  11^,. 
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(4.4)Proposition.  For  any  compactly  supported  ( measurable )  function  <p  on  IRd  and  any 
f  €  tt,  the  following  conditions  are  equivalent: 

(a)  /  e 

(b)  = 

(c) 

Proof.  [BR8x]  use  Poisson’s  summation  formula  (as  did  [SF73],  [DM83],  etc). 
But  since  the  proposition  impose-  no  continuity  requirement  on  p,  the  formula  cannot  be 
applied  directly.  In  this  sense,  the  argument  for  the  corresponding  Proposition  2.2  in  [B87] 
is  incomplete.  This  technical  point  is  handled  in  [BR8x]  by  proving  (b)  to  hold  when  applied 
to  test  functions.  This  makes  the  proof  valid  even  when  p  is  only  a  compactly  supported 
distribution.  Here  is  the  proof,  for  completeness. 

For  any  compactly  supported  test  function  u, 

*'  /)(“)  =  Y  “  a)(“)  =:  Y  M a )’ 

a  a 

with  tp  :  x  >-*  f(x)p(-  -  x)(u)  also  a  compactly  supported  test  function,  and  therefore,  by 
Poisson’s  summation  formula, 


(<P  *'/)(«)  =  ^0(2xa), 

a 

(4.5) 

with 

0  =  f(-iD)(<p(— )u). 

On  the  other  hand, 

(4.6) 

W  */)(«)  =  J  j  <p(v  ~  x)f(x)dxu{y)dy  =  0(0). 

(4.7) 

If  now  /  €  7rv,  then  D0f  €  11^,  for  all  /?,  hence  0(27t a)  =  0  for  all 
(4.6)  and  the  Leibniz-Hormander  identity  (cf.  Appendix), 

a  6  7Ld\ 0  since,  by 

0(0  =  Y  {D0f(-iD))<p{- 0  (-iDfu(0/P'.. 

(4.8) 

Therefore,  from  (4.5)  and  (4.7), 

(<?*'  /)(«)  =  0(0)  =  (¥>*/)(«). 

showing  that  (a)  =>■  (b).  The  implication  (b)  =>  (c)  is  trivial  since  p  *  jr  C  t.  Finally,  if 
(c)  holds,  then  the  linear  functional  2  >— ►  (<^*'  /)(u)  has  support  only  at  the  origin.  Since 
the  collection  of  linear  functionals  {u  >->■  (-iD)0u(2xa)  :  0,a  6  is  globally  linearly 

independent  over  the  space  of  compactly  supported  test  functions,  this  implies  with  (4.5) 
and  (4.8)  that  /  and  all  its  derivatives  must  belong  to  n^,  thus  showing  that  (c)  => 
(a)-  * 

One  recovers  [B87;Prop.2.2]  from  (4.4)Proposition  with  the  aid  of  the  observation 
(already  made  there)  that  p  *’  f  and  /  *'  p  agree  on  TLd  and  f  *'  p  is  a  polynomial  in  case 
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/  €  Jr,  hence  agrees  with  p  *'  f  iff  p  *'  f  £  v.  More  importantly,  (4.4)Proposition  brings  in 
an  important  point  that  simplifies  considerably  the  construction  of  quasiinterpolants,  viz. 
the  equality 

ip*1  =  tp*  on  (4.9) 

which  suggests  the  construction  of  quasiinterpolants  of  the  form 

Qn  :=  p  *'  p*  :  /  *-  Y  ¥>(•  ~  a  W(a  ~  ‘)  (4-10) 


with  the  compactly  supported  distribution  p  chosen  so  that  p*  represents  the  inverse  of  <p* 
on  7rv.  I  pursue  this  point  in  the  next  section. 

As  pointed  out  in  [BR8x],  (4.4)Proposition  has  a  ready  extension  to  exponential  /, 
i.e.,  to 

/  €  Expx  :=  Y,  efl7r 

fl£T 

for  some  finite  T  C  Cd ,  with  eg  :  x  >-*  e9x.  It  is  natural  to  consider  such  exponential  f  i ce 
the  essential  part  of  the  space 


H{p)  :=  {/  €  S(p) :  /  entire  } 
consists  of  exponentials.  I  will  use  the  abbreviation 

o(v>) 

for  the  spectrum  of  II(p),  i.e.,  the  smallest  T  for  which  H(p)  C  Expx.  The  appropriate 
generalization  of  tt^,  is  the  space 


■=  Y 

with  the  largest  f'-invariant  subspace  of 

11^(0)  :=  {p  €  7r  :  p(-iD)p  =  0  on  -  iO  -f  (27rZti\0)}. 

(4.11)Theorem  [BR8x].  Let  p  be  a  compactly  supported  function,  and  let  f  £  Expx. 
Consider  the  following  conditions: 

(a)  /  €  Hv; 

(b)  p*'f  =  p*f- 

(c)  p*'  /  £  Expx. 

Then  (a)  =>  (b)  =>  (c).  If,  in  addition, 

(T  —  T)  n  2ni7La  =  {0},  (4.12) 
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then  (c)  =►  (a)  as  well. 

Here  is  a  proof  outline:  The  implication  (a)=*(b)  follows  from  (4.4)Proposition  by 
shifting  in  the  frequency  domain,  using  the  fact  that  (a)  implies  that  T  C  hence  that 

/  =  fe  with  fg  =  egqe  for  some  qg  6  xv{0).  It  follows  that  each  such  qg  lies  in 

^e_#¥>(0).  Therefore,  by  (4.4)Proposition, 

<P*'  fe  =  eg  {{e-gip)  *'  q6)  =  eg  ((e_*v>)  *  qg)  =  <p  *  fg 
for  each  9  E  0(<p),  and  (b)  follows. 

For  the  implication  (c)  =>■  (a),  decompose  /  into  its  separate  exponential  terms, 
/  =  StfgT  fg,  with  fg  =  egqg.  Assuming  (4.12)  and  (c),  it  is  possible  to  provide,  for  each 
9  E  T,  a  polynomial  p  =  pg  so  that  p(E)f  =  fg  and  p(E)(<p*'  f)  E  egx.  Therefore 

egir  3  p(E)  (<p  *'  f)  =  <p  *'  p(E)f  =  <p*'  fg. 

But  this  says  that  r  :=  e-g(ip  *'  fg)  E  w,  i.e.,  (e_«<p)  *'  qg  =  r  E  it,  therefore  qg  E  ^e_9V(0) 
by  (4.4)Proposition,  i.e.,  fg  E  Hv.  4 

The  assumption  (4.12)  is  essential  here:  If  9,9  E  T  and  0  —  96  2iri2Zd\0,  then 
/  :=  egq  —  e#q  vanishes  on  7Ld  for  any  q  6  tt,  hence  <p  *'  f  =  0,  yet  /  /  0,  hence  does  not 
belong  to  Hv  if  q  has  sufficiently  high  degree. 
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5.  The  construction  of  quasiinterpolants  for  S(<p) 


This  section  follows  closely  the  development  in  [BR8x]. 

In  Section  2,  we  discussed  the  construction  of  F- quasiinterpolants  for  S(ip)  in  the  form 

f*'A  (5.1) 

with  A /  :  x  >-*■  A /(•  +  x )  and  F  an  fJ-invariant,  hence  ^'-invariant  subspace  of  S(y?).  I 
pick  up  on  this  discussion  now  in  order  to  give  a  unified  view  of  the  various  concrete  quasi¬ 
interpolants  already  available  in  the  literature  (cf.,  e.g.,  [SF73],[BH82],  [DM83],  [DM84], 
[DM85],  [DM8x],  [BJ85],  [CJW87],  [CD87],  [CD88],  [CL87],[B87],  [DR8x],  [R8x],  [J88i], 
[J8x],  [BeR8x].  This  is  made  possible  by  the  observation  (proved  in  the  preceding  section) 
that 

ip*'  =  f  *  on  Hy.  (5.2) 

This  suggests  that,  in  trying  to  invert  f*'  on  H^,  we  might  as  well  consider  the  simpler 

problem  of  inverting  the  convolution  <p*  on  H^,,  as  its  inverse  is  given  again  by  a  convolution. 

Thus,  we  look  for  quasiinterpolants  of  the  form 


Qn-=V*'p*  (5.3) 

(instead  of  the  form  (5.1)),  with  p*  any  convenient  convolution  which,  on  F  C  H agrees 
with  T~l.  Here,  as  before,  it  is  assumed  that 

T  :=  (f  *'  j F 


is  1-1,  hence  invertible. 

One  recovers  the  earlier  formulation  (5.1)  for  the  specific  choice  p  \  f  A /(  — •)  (since 

p*f(x)  =  \(E*f)). 

Consider  now  the  problem  suggested  by  the  formulation  (5.3)  in  light  of  (5.2):  For 
given  ^-invariant  F  C  H v  n  S(f>),  find  distributions  p  (of  some  desirable  form)  so  that 
p* |f  =  (vj*|f)-1-  In  fact,  there  is  useful  additional  freedom  here:  It  is  sufficient  to  construct 
p  so  that  p*  =  (0*|f)-1  f°r  some  distribution  ^  for  which  0*  =  f*  on  F. 

For  example,  one  might  choose  to  use  the  particular  distribution 

<P\  :  /•-£  ¥<«)/(«) 

a 

in  place  of  f.  This  is  legitimate:  For,  as  one  computes, 

(<p\  *  /)(*)  =  f\{f{x  -  •))  =  J]  -  a)  =  /  *'  f=f*'  /, 

or 

with  the  last  equality  holding  [B87]  for  any  /  €  S(<p).  Note  that  the  Fourier  transform  of 
f  I  is  the  symbol  of  f,  i.e.,  the  trigonometric  polynomial 

(s?ir= <?’•= 
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which  made  its  appearance  already  in  Section  2. 

We  now  describe  some  concrete  approaches  to  the  construction  of  suitable  p.  While 
the  practical  interest  (if  any)  would  center  on  F  C  it  is  no  complication  to  consider, 
more  generally,  F  =  eg P  C  Hv  fl5(^). 

(i)  Matching  of  Fourier  transform  In  effect,  we  are  looking  for  a  distribution  p  so 

that 

0  *  p  *  /  =  /  -  (5.5) 

for  all  /  in  some  finite-dimensional  exponential  space  F.  By  going  to  Fourier  transforms, 
(5.5)  gives 

=  7 

which  shows  that  we  could  use  here  any  p  for  which  0/2- 1  vanishes  to  sufficiently  high  order 
on  the  (necessarily  finite)  spectrum  of  the  exponential  /.  We  now  make  this  observation 
precise. 

Consider  first  the  case  that  F  is  the  exponential  space  e$P  for  some  F-invariant, 
hence  D-invariant,  polynomial  space  P.  One  computes  (cf.  the  Appendix  on  Polynomial 
Identities)  that 

0  *  (e4D  =  ee  ((0c-«)  *  0°)  =  e«  U°~y[  -  iDfk-iO)  =  M  •  -iD]atj>(-id).  (5.6) 

Given  that  P  is  D-invariant,  this  implies  that,  for  any  polynomial  p  6  P, 
e_g0*(eap)  =  p(--iD)jj(-id)  =  ^Dyp  [  -  tDj70(-ifl)  €  0(-t'0)p  +  P  0  7r<degp.  (5.7) 

7 

In  particular,  0*  maps  eeP  into  itself,  and  is  invertible  on  esP  if  and  only  if  0(-i0)  ^  0. 
Further,  since  p(--iD)  1  =  p,  the  first  equality  in  (5.7)  (with  0  replaced  by  0*p)  shows  that 
0  *  p*  =  l  on  e$P  if  and  only  if  p(-  -  :!/})( 0/2  -  1)(— id)  =  0  for  p  €  P.  This  last  condition 
is  equivalent  to  0/2  -  1  having  a  P(-t-)-fold  zero  at  -id  (i.e.,  p(-z\D)(0/2  -  1)(  —  id)  =  0  for 
all  p  €  P),  since 

p(--iD)  =  J3[r(Z?»(-iD) 

7 

and  P  is  D-invariant.  This  proves  the  following. 

(5.8)Proposition.  If  the  compactly  supported  function  0  satisfies  ifr(-id)  ^  0,  and  the 
finite-dimensional  polynomial  space  P  is  E-invariant,  then  0*  maps  e$P  1-1  onto  itself, 
and  any  convolution  p*  with 

p(-iD)p(-id)  =  p(-i£>)(l/0)(-i0)  V(p  6  P)  (5.9) 

provides  the  inverse  of  0*  on  egP. 

In  applications,  the  polynomial  space  P  is  often  not  known  precisely,  but  its  degree 
can  be  ascertained,  i.e.,  a  k  with  P  C  tt*  can  be  found.  In  that  case,  one  would  satisfy  (5.9) 
for  7 Tfc  rather  than  P,  i.e.,  one  would  make  certain  that  all  derivatives  of  order  <  k  of  p  at 
—id  match  those  of  1/0  there. 
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By  choosing  p  so  that  (5.9)  is  satisfied  with  P  =  P$  for  every  0  G  0(</j),  one  obtains 
a  suitable  distribution  p  and  a  quasi-interpolant  p  *'  p*.  For  example,  if  p*  =  q(iD)  for 
some  polynomial  q,  then  fi  =  q,  while  the  choice  p*  =  q{E)  leads  to  the  ‘trigonometric’ 
polynomial  fi  =  g(e‘^).  In  either  case,  appropriate  osculatory  interpolation  to  l/p  at 
—iQ(p)  provides  an  appropriate  p.  In  the  first  case,  p  is  a  linear  combination  of  values 
and  derivatives  at  the  origin,  while,  in  the  second  case,  p  employs  only  function  values  at 
some  points  from  7Z3.  More  generally,  one  could  use  p  of  the  form  /  Erex  Px{iD)f(x) 
in  which  the  polynomial  px  is  to  be  chosen  so  that  the  exponential  fi  =  ErexP*()e,a^' 
osculates  to  l/p  appropriately  at  —  i@(<p). 

If  one  uses  p  of  the  form  u>|,  then  p*  =  in(a)[a]  is  a  difference  operator,  hence 
commutes  with  <p*'  and  thus 


V  *' (m  *  /)  =  (m  *  V3)  *'  /• 

This  provides  us  with  a  quasi-interpolant  of  the  simple  form  ip*',  with  ip  €  S(<p),  and  with 
the  support  of  ip  not  exceeding  the  sum  of  the  supports  of  p  and  <p.  In  fact,  it  is  shown  in 
[R8x]  that  the  support  of  the  difference  operator  p*  can  be  chosen  so  that 

diam  supp  ip  <2  diam  supp  <p, 

in  contrast  to  the  minimal  polynomial  procedure  below  in  which  the  inverting  difference 
operator  is  supported  on  a  relatively  large  domain. 

We  summarize  the  Fourier  transform  approach  in  the  following 

(5.10)  Theorem.  Let  p  be  a  compactly  supported  function  whose  Fourier  transform  does 
not  vanish  on  let  F  be  an  E-invariant  subspace  of  fl  S(p)  and  p  a  compactly- 

supported  distribution.  Then 

Qn  :=  V3  *'  p* 

is  an  F-quasiinterpolant  (i.e.,  is  the  identity  on  F)  if  and  only  if 

p(-iD)fi(-iO)  =  p(-iD)(l/ip)(-i0),  alleepe  F, 

where  ip  is  any  compactly  supported  distribution  for  which  ip*  coincides  with  ip*'  on  F. 
Suitable  choices  for  ip  are  ip  =  <p  and  ip  =  <p\. 

By  comparison  with  the  Fourier  transform  approach,  the  remaining  two  approaches 
don’t  seem  very  efficient.  They  apply  directly  only  to  F  of  the  simple  from  egP. 

(ii)  Minimal  polynomial  Here,  one  would  choose  an  ‘easily  computable’  distribution 
ip  for  which  ip*  =  <p*  on  F  and  observe  that,  since  T  —  {<p*')\f  is  invertible,  we  can  represent 
T_1  as  p(ip*)  for  some  polynomial  p  (i.e.,  obtain  the  inverse  of  the  operator  ip*  as  a  linear 
combination  of  powers  of  ip*).  We  obtain  such  a  polynomial  as  a  multiple  of 

(mT(0)  -  %)/()*, 

with  mr  the  minimal  (annihilating)  polynomial  for  T.  It  may  be  hard,  in  general,  to 
produce  this  polynomial,  particularly  if  the  space  H ^  0  5(i^)  is  not  known  precisely.  But, 
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we  conclude  from  (5.7)  that,  on  the  exponential  space  cjtt,  (p *  is  degree- preserving  in 
the  sense  that,  for  any  polynomial  p, 

f  *  eeP  €  <p{-i0)(egp )  +  e97r<degp. 


In  fact,  (5.7)  implies  that 


<P*eepe  <p(—iO)(egp)  +  e«7rdegp_fc 

in  case  <p  —  <p(-id)  has  a  zero  of  order  k  at  —id.  For  example,  if  cp  is  radially  symmetric, 
i.e.,  <p(- 1)  =  <p(x),  then  <p(x)  =  <p(0)  +  0(|xj2),  hence 

<P  *  p  €  £(0)p  +  7Tdegp_2  (5.11) 

in  that  case,  as  was  already  pointed  out  in  [CD87].  In  any  case,  this  makes  available  the 
linear  polynomial  rs  :=  ip(-id)  -  •,  for  which  rg(<p*)  is  degree-reducing  on  e$7r.  It  follows 
that,  for  any  F  in  egir,  there  is  a  suitable  power  (r$)n  of  r  which  annihilates  T,  yet  does 
not  vanish  at  0  (since  rg(0)  =  <p(— id)  ^  0  as  T  is  assumed  to  be  invertible)  and  therefore 
is  available  for  the  construction  of  T-1  in  the  form  p(ip*). 

(iii)  Recurrence  Equation  (5.7)  suggests  the  solution  of  the  equation  ip*l  =  /  €  e$P 
by  backsubstitution,  i.e.,  by  recurrence,  since  it  implies  that,  for  /  =  e#p  €  e#P, 

therefore  (using  the  invertibility  of  ip*  on  egP ) 

/  =  -  wr^(-id), 

7 

hence 


(0*)  1/  =  (V>*)  l(e*p)  =  (/-J>*) 

7/0 


-l 


(e9F>7p)[  -  iD]1  i>{-id)  /ip{-id).  (5.12) 


For  a  general  exponential  /,  the  resulting  solution  depends  of  course  on  the  choice  of 
ip*,  but  necessarily,  since  ip*  =  <p*'  on  F,  this  solution  is  independent  of  ip  in  case  /  €  F 
(as  it  then  equals  (<p*')_I /).  Note  that  the  recurrence  does  require  the  numbers 

D^’p(-id) 


for  all  7  for  which  F7p  ^  0.  It  also  requires  the  generation  of  the  ‘tree’  of  derivatives  of  p 
which,  during  the  solution  process,  would  have  to  be  traversed  in  the  order  opposite  to  the 
order  of  its  generation. 

By  contrast,  the  use  of  Appell  polynomials,  detailed  in  Section  2,  does  not  require 
explicit  information  about  the  derivatives  of  ip.  In  effect,  it  computes  the  solution  pa  for 
the  equation  ip *?  =  []°  for  all  a  needed  directly  (by  simple  recurrence),  and  obtains  the 
solution  for  general  p  in  the  form  ^2apaDap( 0). 
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6.  Projectors  from  quasiinterpolants 


Linear  projectors  R  onto  5  are  particularly  desirable  approximation  maps  into  5  since 
they  provide  the  estimate 

|A/  -  XRf\  <  dist(A,  ran  i?,)||l  -  i?||  dist(/,  5) 

for  any  function  /  and  any  linear  functional  A.  (Here,  ran  R'  consists  of  the  interpolation 
conditions  for  R,  i.e.,  raniZ'  =  {A  :  A/  =  A/Z/,  V/}.)  In  particular,  if  we  are  able  to 
construct  a  bounded  sequence  {Rk)  of  projectors  for  the  directed  set  (5a),  then  we  know 
that 

dist(/,  Sh)  =  0{\\f  -  Rhf\\), 

hence  can  determine  the  approximation  order  of  (5a)  from  the  behavior  of  || /  -  ZZa/||  for 
smooth  /  as  |/i|  — ►  0.  Such  behavior  will  be  hard  to  determine  unless  R a  is  local.  Thus  we 
are  looking  for  good  quasiinterpolants  which  are  also  projectors. 

Note  that  the  existence  of  such  a  good  projector  implies  that  the  approximation  order 
is  at  least  m  in  case  7r<m  C  5a  for  all  h.  Thus  we  cannot  expect  to  construct  such  good 
projectors  for  spaces  like  Cl  cubics  on  the  three-direction  mesh,  not  even  for  its  locally 
spanned  subspace  5({A/i22i  M212})  (cf.  Section  3). 

Neither  can  we  expect  to  construct  such  a  good  projector  for  5  unless  we  have  available 
a  basis  $  for  5  which  is  locally  finite  (i.e.,  locally  finite-dimensional)  in  the  sense  that, 
for  all  n,  at  most  finitely  many  <p  €  $  have  some  support  in  any  particular  open  ball.  With 
Bn  :=  Bn( 0)  the  ball  of  radius  n  around  the  origin,  this  is  equivalent  to  saying  that,  for 
every  n,  the  set 

t'n  :=  {<p  €  $  :  supp  <p  Pi  Bn  ^  0} 

is  finite.  For  any  such  collection,  it  makes  sense  to  define  the  linear  map 

$  :  IB*  — ►  F  :  c  $c  :=  ^  p>  c(p>), 

with  the  (possibly)  infinite  sum  taken  pointwise,  or,  more  strongly,  uniformly  (since  triv¬ 
ially)  convergent  on  compact  sets.  (There  seems  to  be  no  harm  in  using  the  letter  $  both 
for  the  set  and  the  map  induced  by  it,  much  as  we  might  use  the  same  letter  A  to  denote 
a  matrix  and  the  sequence  A  =  [al,...,an]  of  its  columns.)  We  say  that  $  is  (globally) 
linearly  independent  if  the  map  $  is  1-1.  In  that  case,  $  provides  a  basis  for  its  range 

5  :=  ran  $. 

This  implies  that  any  linear  man  into  5  can  be  written  in  the  form 

*>€♦ 

for  suitable  linear  functionals  Av  (and  with  v?A  :  /  <pA(/)).  In  particular,  such  a  map  is 

a  linear  projector  if  and  only  if,  for  every  ip  €  A^,  extends  the  coordinate  functional 

A”5 -*  IR  :  /  ~  (*-7)(¥>). 
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Consequently,  such  a  projector  cannot  be  local  unless  each  coordinate  functional  is  local, 
i.e.,  satisfies 

f\A  =  0  =»  A  If  =  0 

for  all  /  £  5  and  some  bounded  set  A.  This  makes  the  following  recent  result  [BeRSx]  of 
Ben-Artzi  and  Ron  particularly  interesting. 

(6.1)Theorem.  Let  $  be  a  locally  finite,  globally  linearly  independent  collection.  Then, 
for  every  there  exists  a  ball  B  -  so  that  ($c) \b  =  0  implies  c(<p)  =  0. 

Proof.  Without  loss,  we  assume  that  (after  a  suitable  translation)  <p  has  the  origin 
in  its  support.  In  view  of  the  definition  of  $c  :=  <pc(<p)  as  a  pointwise  sum, 

ker  $  =  Q  ker  Rn , 

n 


with 

Rn  :  IR*  —» IRS"  :  c  •-»  <pc(<p))\Bn- 

We  claim  that  the  sequence  ker  R\  D  ker  R2  D  •  •  •  is  decreasing  fast  enough  so  that,  for 
some  n ,  [ip]  =  0  on  kerf?n,  i.e.  c  £  keriZn  implies  c(<p)  =  0.  The  proof  is  by  contradiction, 
i.e..  by  showing  that,  in  the  contrary  case,  ker$  =  f)n  ker  Rn  ^  0. 

Whether  or  not  the  sequence  c  belongs  to  ker  Rn  depends  entirely  on  its  behavior  on 
vn,  i.e.,  c  £  ker  Rn  iff  rnc  £  ker  Rn,  with 

rn  :  c  •—  Xvnc 

the  natural  projection  of  IR*  onto  the  space  of  coefficient  'sequences'  having  support  only 
on  i/n.  Thus,  with 

vo  '■= 

which  is  consistent  with  the  otherwise  evident  relationship 

Vj\  C  Vrc, 

since  <^(0)  /  0  by  assumption,  we  can  describe  the  contrary  case  as  asserting  that,  for  every 
n,  r0  is  not  trivial  on 

A/„  :=  ker  Rn  ft  ran  rn, 

and  this  linear  space  is  finite-dirncnsionnlby  the  local  finiteness  of  $.  Since  rn(.\/,l+1 )  c  Mn 
for  any  n,  this  implies  that,  for  any  n  and  k, 

A fn,k  ■—  j  ( ) 

is  a  nontrivial  subspace  of  A/n  since,  by  assumption,  ker  Rn+k  contains  some  c  with  c(<^)  ^  0, 
hence  c'  :=  rn  ■  ■  •  rn+*c  £  Mn\ 0  since  c'(p)  =  c(p)  ^  0.  It  follows  that  Mn  0.Mn, i,  --  is  a 
decreasing  sequence  of  nontrivial  finite-dimensional  spaces,  hence  has  a  nontrivial  limit 

f^n  ■ ~  P|  A/n>fc  =  Pi  r„A/n+ijt_i  —  rnI\n+l. 
k  k>  1 
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Therefore,  starting  with  some  Co  G  /<o\0,  we  can  find,  by  induction,  the  sequence  co,Ci,C2, 
...so  that  cn  €  Kn  and  rn_i cn  =  cn_j  for  n  =  0,1,2, —  Its  limit,  c^,  provides  a 
nontrivial  element  of  ker$  =  finkeri2„,  thus  contradicting  the  linear  independence  of  $. 

♦ 

We  call  A  for  which  f\ a  —  0  =>  Xf  =  0  for  all  /  G  5  a  determining  set  for 
A  G  S'.  Note  that  a  linear  functional  on  a  function  space  on  IR^  may  have  many  different 
determining  sets.  For  example,  the  coordinate  functional  of  a  B-spline  has  every  open 
subset  of  the  B-spline’s  support  as  a  determining  set. 

The  theorem  implies  that  each  A°  is  local,  hence,  by  the  local  finitenes  of  5,  sees  5  as 
a  finite-dimensional  linear  space,  i.e.,  the  space 

Sa  :=  SjA 

for  some  bounded  A.  This  implies  that  each  A°  is  bounded  (in  whatever  norm  we  care  to 
impose  on  5),  hence  admits  extensions  to  a  continuous  linear  functional  with  support  in  A 
on  whatever  normed  linear  superspace  X  we  wish  to  approximate  from  5(9). 

What  is  offhand  missing  is  the  uniform  boundedness,  i.e.,  the  finiteness  of 

sup  !|A°||/!MI, 

vS* 

or,  equivalently,  the  finite  condition  of  the  basis  $.  Such  information  is  needed  for  the 
construction  of  bounded  linear  projectors  to  5.  We  do  obtain  such  boundedness  when 
$  =  {<£(•  -  a)  :  a  6  7Ld}  and,  correspondingly, 

5  =  5(9) 

for  some  compactly  supported  function  9.  For,  in  that  case,  the  fact  that  Eaip*c  -  9*(  Eac) 
implies  that 

\0  _  \0  pa 

Av(-a)  ~  V  > 

hence  that 

/  =  v3*'A0*/  V(/G5(9)) 

with 

A0:/~A»/(~) 

bounded.  Each  Qs  '■=  9  *'  A *  so  obtained  is  a  uniformly  local  linear  projector  onto  5(9). 

Consider  now  the  actual  construction  of  such  A,  in  the  spirit  of  [J Sx] ,  [J88],  and 
[BeR8x],  and  in  imitation  of  the  dual  functionals  for  the  univariate  B-splines  discussed  in 
the  first  section.  This  requires  us  to  consider  the  linear  functionals  of  the  form 

[x]p(D):f~(p(D)f)(x) 

for  x  6  IR1*  and  p  €  w.  We  will  write 


p*  :=  (0 }p(D)\Sa- 

Specifically,  assume  that,  for  some  1  G  A  and  some  .E-invariant  polynomial  space  P.  the 
map 

P-~(SAY:p>~  [x]p(D)  (6.2) 
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is  invertible.  Then  there  exists  a  unique  p  €  P  so  that 

Ao  =  [x]p(D)  on  SA. 

If  now  also  F  is  a  D-invariant  (hence  smooth)  space  for  which 

F\a  -  SA  =  S(tp)\At 
then,  for  any  y  6  lRd,  F  =  EyF,  hence 


P->F':p~[y]p(D)  (6.3) 

is  invertible  for  each  y  €  IRd,  thus  providing  a  unique  extension  of  Ao  of  the  form 

Ao  =  [x]Px(D)  on  F 

for  every  x  €  IR*.  In  order  to  discover  the  behavior  of  px  as  a  function  oi  x,  take  any  basis 
B  for  F,  let  ( Pb*)beB  be  the  corresponding  dual  basis  in  [0]P(Z>),  and  let  p  be  an  arbitrary 
linear  functional  on  SA.  Then 

P  =  X>P6*  =: 

b6B 


therefore 

M  =  0 uE~x)E*  =  fK  ~  *)  [*]?*(£), 

b 

using  the  fact  that,  for  any  x  and  any  smooth  /  and  any  p  €  tt, 


p*(Exf)  =  [0)p(  £>)/(•  +  x)  =  [x]p(D)f. 

This  extracts  the  essential  features  of  the  argument  in  [BeR8x]  for  the  following  proposition, 
which  generalizes  the  dual  functional  formula  (1.6)  for  univariate  B-splines. 

(6.4)  Proposition.  If  $  is  globally  linearly  independent  and  F  is  a  D -invariant  space  which 
agrees  ■■■ <’ith  ran4>  on  some  determining  set  A  for  A0  :=  A°  for  some  y>  G  then,  for  an 
arbitrary  basis  B  of  F  and  for  every  x  €  IR^,  the  linear  functional 


J2<l(D)b(-x)lx}Pb(D)  (6-5) 

beB 

extends  Ao|f-,  where  ( pb )  is  the  unique  ‘dual  basis’  for  B  in  any  E-invariant  polynomial 
space  P  for  which  the  map 

P  —  F'  :py^  [0]p(£>)|F 

is  invertible,  and  q  is  the  unique  element  of  P  for  which  [0]qr(Z?)  =  Aq. 


Since  (6.5)  only  extends  Ao|f,  it  extends  Ao  only  if  x  is  in  some  determining  set  A  for 
A0  on  which  F  =  ran$.  Thus  we  could  recover  the  formula  (1.6)  from  this  proposition 
(using  F  =  P  =  7T<*.).  On  the  other  hand,  the  fact  that  any  nontrivial  knot  interval  in  the 
support  of  a  univariate  B-spline  is  determining  for  its  coordinate  functional  is  most  easily 
proved  with  the  aid  of  the  dual  functional  formula  (1.6). 

The  propositon  makes  no  use  of  any  local  linear  independence.  (To  recall,  $  is  locally 
linearly  independent  if  {ip \A  :  p>  €  $,<P\A  ^  0}  ’s  linearly  independent  for  any  open  set 
A.)  While  it  is  an  open  question  whether  the  linear  independence  of  $  implies  the  linear 
independence  of  $n  for  all  sufficiently  large  n  (the  guess  is  that  it  does  not),  there  are  simple 
examples  to  show  that  linear  independence  does  not  imply  local  linear  independence.  For 
example,  the  integer  translates  for  the  characteristic  function  of  the  interval  [0,1/2]  are 
linearly  independent  but  not  locally  linearly  independent. 
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7.  Good  projectors  into  smooth  bivariate  polynomial  spaces 

Once  we  consider  spline  spaces  more  general  than  those  generated  by  the  translates  of 
a  single  compactly  supported  function,  we  are  at  the  frontier  of  research.  The  only  case 
studied  in  any  depth  is  that  of 

S  =  <A 

of  all  functions  in  Cp  which  are  piecewise  polynomial  of  degree  k  at  most,  with  A  a  trian¬ 
gulation  of  some  domain  in  IR2.  The  results  are  as  follows. 

It  has  been  known  for  some  time  (see,  e.g.,  [Z74]  and  references  there)  that  5  has  full 
approximation  order,  i.e., 

dist(/,S)  =  0(|A|*+1) 

for  all  smooth  functions,  in  case  k  >  4p  +  1,  since  it  is  then  possible  to  construct  a  good 
projector  by  local  Hermite  interpolation  in  each  triangle.  Here,  jA|  is  the  largest  diameter 
of  a  triangle  in  A.  (I  use  here  the  adjective  ‘full’  since  the  approximation  order  from  the 
entirely  unconstrained  (and  larger)  space  tt^a  cannot  be  better  than  0(|A|*),  by  arguments 
like  those  we  used  in  the  univariate  case  in  the  first  section.) 

It  was  shown  in  [BH88]  that  full  approximation  power  was  had  by  5  even  if  only 
k  >  3p  +  2.  The  argument  given  is  not  constructive,  but,  as  is  made  clear  in  [B90], 
contains  the  essential  information  for  the  construction  of  a  stable,  local  basis,  hence  of  a 
good  projector,  for  S.  [CL88]  provide  an  explicit  construction  of  a  good  projector  onto 
a  subspace  of  5  which  contains  hence  already  achieves  the  full  approximation  order. 
These  results  are  sharp,  i.e.,  already  for  A  the  three-direction  mesh  (cf.  [BH88]  for  p  <  3, 
R.-q.  Jia  for  general  p),  the  approximation  order  is  less  than  k  +  1. 

Yet  these  results  have  a  puzzling  aspect.  While  the  local  basis  constructed  consists 
of  ‘vertex  splines’,  i.e.,  of  splines  with  support  within  the  triangles  sharing  one  vertex,  the 
stability  of  the  basis,  i.e.,  its  condition  number,  depends  on  the  mesh  A.  This  is  to  be 
expected  since  there  should  be  difficulties  with  skinny  triangles,  i.e.,  triangles  with  a  small 
angle.  What  is  annoying  is  that  there  are  also  difficulties  when  A  has  many  near-singular 
vertices,  i.e.,  vertices  with  just  four  edges  and  with  opposing  edges  nearly  parallel.  Jia  has 
recently  established  the  existence  of  a  local  basis  for  a  for  k  >  3p  +  2  whose  condition  is 
bounded  in  terms  of  the  smallest  angle  in  the  partition  only,  but  some  of  its  elements  may 
fail  to  be  vertex  splines. 

A  good  and  challenging  test  case  is  provided  by  the  space  of  bivariate  C^-cubics.  For 
A  the  three-direction  mesh,  its  approximation  order  is  known  to  be  3  rather  than  4  (cf. 
Section  3),  but  it  is  not  known  what  its  generic  approximation  order  is.  Since  there  are 
partitions  A  in  which  some  triangles  are  not  in  the  support  of  any  compactly  supported 
element,  hence  for  which  5  fails  to  contain  a  local  partition  of  unity,  my  guess  is  that, 
generically,  bivariate  C*-cubics  have  approximation  order  0.  In  this  regard,  my  working 
hypotheses,  in  order  of  increasing  content,  are  the  following.  While  they  are  meant  for  any 
‘directed’  set  (5*)  with 

<,a„ 

and  |h|  :=  |A/,|,  one  would  try  them  first  for  scales,  i.e.,  for  A^  =  hA  for  some  fixed  trian¬ 
gulation  A  for  IR2.  Also,  local  support  means  ‘support  of  0(|A/,|)\  and,  correspondingly, 
local  partition  of  unity  or  local  basis  means  a  partition  of  unity  or  a  basis  whose  elements 
have  local  support. 
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Conjecture  1.  If  (Sh)  has  positive  approximation  order,  then  5/,  contains  elements  with 
local  support. 

Conjecture  2.  If  the  scale  (Sh)  has  positive  approximation  order,  then  Sh  contains  a  local 
partition  of  unity. 

Conjecture  3.  The  approximation  order  of(Sh)  equals  that  of  (S*'0'),  with 

Sloc  :=  span{/  €  S  :  supp  /  compact). 

Conjecture  4.  The  approximation  order  of  (Sh)  can  always  be  realized  by  a  corresponding 
set  (Qh)  of  good  quasiinterpolants  ( with  ran  Qh  C  Shoc  necessarily). 

A  first  stab  at  Conjecture  2  was  made  in  [BD85],  in  a  univariate  context,  with  S 
restricted  to  be  ^-invariant.  [J882]  extended  this  to  a  univariate  result  which,  if  true  also 
in  the  multivariate  setting,  would  give  the  following  positive  answer  to  Conjecture  4. 

Conjecture  5.  If  S  is  an  E-invariant,  locally  finite-dimensional  space  of  functions  closed 
under  uniform  convergence  on  compact  sets,  then  the  corresponding  scale  (Sh)  has  approx¬ 
imation  order  m  if  and  only  if  S  contains  a  compactly  supported  function  if  whose  Fourier 
transform  satisfies  the  Strang-Fix  conditions  of  order  m,  i.e.,  ^>(0)  =  1  and  Datf  =  0  on 
2x7Zd\0  for  all  |o|  <  m. 


Free  knot  splines 

In  these  lectures,  I  have  failed  entirely  to  cover  the  multivariate  equivalent  of  splines 
with  free  knots.  This  is  already  a  challenging  problem  in  the  tensor  product  case,  which, 
however,  doesn’t  strike  me  as  germane  since  the  meshes  there  have  global  effect  while  an 
adaptive  meshing  should  have  local  adaption  as  its  goal.  There  is  an  early  paper  [BRi82] 
on  the  local,  adaptive  refinement  of  a  rectangular  grid.  While  [BIti82]  deals  only  with 
piecewise  polynomials  without  smoothness  constraints,  a  corresponding  development  for 
smooth  bivariate  piecewise  polynomials  spanned  by  simplex  splines  can  be  found  in  [D82]. 
Both  of  these  papers  (and  others)  deal  with  adaptive  refinement.  The  problem  of  the 
construction  of  an  optimal  partition  A,  i.e.,  a  partition  of  the  given  domain  G  into  a  given 
number  of  triangles,  say,  so  as  to  minimize  dist(/,  7t£a)  is  considerably  harder.  There  is 
Nadler’s  result  [N85]  on  the  locally  optimal  shape  of  such  triangles  when  k  =  1  and  p  =  0. 
A  totally  different,  and  very  intriguing,  approach  has  been  started  by  DeVore  and  Popov 
[DP87]:  They  determine  a  best  approximation  (in  the  Lp-norm)  from  the  space  of  linear 
combinations  of  n  characteristic  functions  of  shifted  (hyper)cubes  of  arbitrary  size.  In 
particular,  they  do  not  insist  that  these  cubes  not  overlap  and  thus  avoid  a  major  difficulty 
associated  with  the  construction  of  optimal  meshes.  Subsequent  work  by  DeVore,  Jawerth 
and  Popov  (in  progress)  makes  use  of  the  existing  and  developing  theory  of  multiscale 
expansions  (wavelets)  to  extend  these  results  to  box-splines  and  other  smooth  piecewise 
polynomials.  Since  this  uses  dyadic  subdivision,  it  does  look  more  like  adaptive  refinement 
than  optimal  mesh  generation. 
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8.  Appendix:  Appell  polynomials 

Traditionally,  an  Appell  sequence  (pa)  of  polynomials  is  characterized  by  the  fact  that 
Po  =  1  and  D^pa  =  pa-p. 

Thorne  [T45]  pointed  out  the  much  more  useful  characterization  in  terms  of  a  lin¬ 
ear  functional  A,  normalized  so  that  A()°  =  1,  with  which  there  is  associated  uniquely  a 
polynomial  sequence  ( pQ )  by  the  conditions  that 

(i)  Pa,  e  TTa  (ii)  (A Dp)pa  =  Sp_a.  (8.1) 

Actually,  Thorne  and  others  only  consider  univariate  polynomials,  but  the  multivari¬ 
able  extension  is  clear  in  view  of  the  now  standard  multivariable  index  notation. 

Standard  examples  include 

A  =  [0],  giving  the  normalized  powers  or  Taylor  polynomials  pa  =  |]“; 

A  :  /  !-►  JqI/Jq  1 ,  giving  the  Bernoulli  polynomials  for  SI. 

A  =  ([0]  +  [l])/2  on  C(IR),  giving  the  Euler  polynomials. 

Existence  and  uniqueness  of  the  Appell  sequence  for  given  A  follows  directly  from  the 
definition.  In  fact  (as  Amos  Ron  kindly  pointed  out),  (ii)  implies  (i):  If  7  is  a  maximal  index 
with  the  property  that  £)'1'pa( 0)  ^  0,  then  Dypa  is  a  nonzero  constant,  hence  A D^pa  ^  0, 
which  implies  that  7  =  0.  With  that, 

=  E  nvr), 

T ,<a 


and  the  condition  A D^pa  =  &a-p  >s  equivalent  to  the  linear  system 

£  W^Mt)  =  b-c  V(/?  <  a), 

■y<a 


which  has  a  triangular  coefficient  matrix  with  unit  diagonal  (since  Af]0  =  1  by  assumption), 
hence  has  exactly  one  solution,  and  this  solution  can  be  obtained  by  backsubstitution, 

p*  =  r  -  E  p-t  Mf-  (8.2) 

0?a 

The  correctness  of  this  formula  can  be  verified  directly  by  induction  on  a: 

A Wpa  =  AD70°r  -  A[f 

=  A[f-7  -  Alf7  =  0 

for  7  <  a,  while  A Dapa  —  AD0}]0  =  A()°  =  1. 

The  uniqueness  implies  that 

DPPa  =  Pa-0 ,  (8.3) 

hence 

p.  =  E  d”~M7) 

7<0r 

for  some  (as  yet  mysterious)  sequence  c. 
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(8.5) 


(8.4)Proposition.  For  <j>  :=  <p\  :/»-»■  A/(--)  and  for  all  a, 

p*  =  l  -iD]°(  i/JjXO)  =  n”'7  [o][  -  .-crt  i/?»). 

*y<a 

Proof.  Convolution  maps  polynomials  to  polynomials.  Precisely,  for  any  p  €  7 r 
and  any  <£, 

<t>*p=p(.~  w)m  =  22DyP[-  iDYm- 

Consequently, 

4>  *  P  e  4(0)p  +  TTcdegp.  (8.6) 

This  implies  that  <f>*  maps  each  wa  into  itself,  and  is  invertible  on  each  such  tt a  iff  <£(0)  ^ 
0.  Moreover,  on  ttq,  the  action  of  <f>*  is  entirely  determined  by  the  numbers  Dp(p{ 0),  (3  <  a. 
In  particular,  if  p  is  any  function(al),  then  convolution  with  the  composite  tp  :=  <f>  *  p  also 
carries  each  ira  into  itself;  it  is  the  identity  on  7ra  iff  the  Fourier  transform  tp  —  <pp  satisfies 

D^(0)  =  S 7  V(7  <  a).  (8.7) 

In  particular,  take  now  convolution  with  any  functional  p  for  which  <p\p  —  1  a-fold  at 


0,  i.e.,  for  which 

D0$xp)(O)  =  60  V(/3  <  a). 

(8.8) 

Then 

D0m  =  D0( l/&)(0)  V(/3  <  a) 

(8.9) 

and  p *  represents  the  inverse  of  <f>\*  on  xa.  In  particular, 

II0  =  <t>\  *  (p  *  0^)  V(/3  <  a). 

(8.10) 

Consequently  (using  the  fact  that  D  commutes  with  convolution), 

Sa-p  =  [O]D0Ua  =  [0]4>X  *  (D0(p  *  0°)), 
while  [0](<^a  *  /)  =  A/.  Since  p  *  []“  €  irai  it  follows  that 

pc = m  *  n°  =  E  fJ'i  -  iur(i/?A)(o). 

y 


* 

The  proof  suggests  that  it  would  be  particularly  helpful  to  identify  linear  functionals 
p  satisfying  (8.8)  since  that  would  obviate  use  of  the  recurrence  for  the  construction  of  the 
Pp  for  (3  <  a.  This  is  the  same  question  as  finding  a  formula  for  the  linear  functional  for 
which 

9a  :=  4>\  *  0“ 

are  the  Appell  polynomials.  That  these  qa  form  an  Appell  sequence  in  the  traditional 
sense  of  the  word  follows  from  the  fact  that,  on  7r,  <p+  commutes  with  differentiation.  This 
suggests  an  alternative  (and  simpler)  way  for  generating  an  Appell  sequence  from  a  linear 
functional.  It  would  be  harder,  though,  to  know  the  interpolation  conditions  which  go  with 
the  expansion  in  terms  of  the  qa. 
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9.  Appendix:  Polynomial  identities 


(d.l)Definition. 

[]n  :  IR  -h.  IR  :  x  >-+  xn/n\ 

[1“  :  IRd  -*  IR. :  i«n[x(j)]“u) 

i= i 

(9.2) Multinomial  identity. 

[*+*+■••+*]"=  £  I'ftr-w' 

i+v+--+<=a 

Proof.  by  induction  on  |o|.  4 

(9.3) Taylor.  For  any  polynomial  p, 

p(x+v)  =  ]ri*r£ap(y)- 

a 

Proof.  For  the  particular  polynomial  p  =  I]*5, 

p(*  +  y)  =  £  [*]>]'’■“  =  £  W“(0“p)(p)- 

a  a 

♦ 

(9.4) Leibniz.  For  any  functions  f,g,. .  .,h  and  any  scalar  s, 

[sD]a(fg  £  \*DYf  [sDfg  ■  ■  •  [sF^/i 

V’+T+-”+7)=o 

Proof.  by  induction  on  |a|.  4 

(9.5) Leibniz-Hormander.  For  any  polynomial  p,  scalar  s,  and  functions  f,g, 

P(»J >)(/»)  =  £  J>)  (»c)/)  Is£)lV 

P 

Proof. 

p(aF>)(/<7)  =  5^a(a)ls£>]“(/5) 

Or 

=  £<‘(“)£[*0]”"'7[»1>]‘’»  =  £  ((|r|"p)  (»D)f)  [nfllV 

a  p  P 

* 

Note.  The  identity  is  linear  in  p,f,g,  hence  verifiable  by  checking  it  just  for  pure 
powers. 
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(9.6)  Convolution.  For  any  compactly  supported  <f>  and  any  p  G  n, 

<t>*P=p(-~  w)m  =  y  irDM-iD)m  =  YD1p  i  -  iDvm- 


Proof. 


hence 


So, 


Dp4>=  J  4>{y){-iy)0 e  xy()dy, 

(-i-D)^(O)  =  J  <t>(y){-y)pdy. 

**na  =  J <t>{- -  y)h}a dv  =  J  <t>{y)l- -y\a dv 

=  [  m  Y  na’p[  -  yfdv  =  E  /  Ml  -  yfdy . 

j  n  a  J 


p 

=  E  -  iDfm  =  i  •  -mam- 

p 
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